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pg ' The He abundance is not constrained by helioseismic data. Mixing of He 

CS ' inside the solar core by processes not included in the solar standard model 

(SSM) has been recently proposed as possible solution to the solar neutrino 
0> , problem. We have performed a model independent analysis of solar neutrino 

^^ I fluxes using practically arbitrary ^He mixing. In addition, we have been si- 

i-^ ' multaneously varying within very wide ranges the temperature in the neutrino 

production zone and the astrophysical factors, Sn and 534, of thep + ^Be and 
^He + ^He cross sections. Seismic data are used as constraints, but the solar- 



c/2 , luminosity constraint is not imposed. It is demonstrated that even allowing 

' ^He abundances higher by factors up to 16 than in the SSM, temperatures 5% 

(or more) lower, the astrophysical factor 5i7 up to 40% higher, and varying 

/\^ • S34 in the range (-20% ,+40%), the best fit is still more than 5a away from the 

j^ ■ observed fluxes. We conclude that practically arbitrary '^He mixing combined 

with independent variations of temperature and cross-sections cannot explain 

the observed solar neutrino fluxes. 



I. INTRODUCTION 

Detected neutrino fluxes in all five solar neutrino experiments (Homestake 0], SAGE 0], 
GALLEX p|, Kamiokande [^], and Superkamiokande 0) are smaller than those predicted 
by SSM's. The status of the solar neutrino measurements is that of disappearance oscillation 
experiments. 

How reliable is this conclusion? In other words, how reliable are the SSM predictions of 
solar neutrino fluxes? 

1. Helioseismic observations of density and sound-speed-squared profiles are in agree- 
ment with SSM predictions throughout the sun with an accuracy better than a fraction 

1 



of percent 0. This agreement was recently found to be valid also in the inner solar core, 
R < O.IRq, where neutrino fluxes are mostly produced JTHTTl. 



In the literature the relevance of this agreement for the neutrino fluxes has been ques- 
tioned with two kinds of objections. 

First, seismic measurements give only the sound speed squared, cj., but not directly the 
temperature T, which is mostly relevant for neutrino production. The connection between 
the uncertainties of c^ and T is 6cl/cj. = ST/T — Sfi/fi, where /x and 6fi are the molecular 
weight and its uncertainty. In principle, one can imagine that 6T/T much larger than Sc1/c1 
could be compensated by correlated changes of the molecular weight. Note, however, that 
this compensation (fine tuning) should work for all distances, and this is harder to imagine. 

Second, seismic uncertainties grow for R < 0.05Rq. Nevertheless, these uncertainties are 
too small to solve solar-neutrino problem. Careful analyses of the inversion method and of 
the available helioseismic data produce c^ in the agreement with the SSM value at level of 
1% at -R = 0, and better in the neutrino production zone |T^ . 

2. Nuclear cross-sections for neutrino production could be responsible for the observed 
neutrino deficit. 

Recent measurements [|1^ of the '^He -|- ^He cross section at the energy of the Gamow 



peak in the sun have eliminated a dangerous source of potential uncertainty. Moreover, even 
before this experiment, it had been already demonstrated that any combination of unknown 
cross sections and variations of central temperature cannot not explain the combined results 



of solar-neutrino experiments [14 



As everything in the world, SSM is not perfect. It does not include, for example, rotation 
and hydrodynamical processes. A bump in the sound-speed profile at i? ~ 0.7 Rq shows sta- 
tistically significant disagreement between SSM and helioseismic observations. Nevertheless, 
this disagreement is far too small to affect neutrino fluxes. It is also most likely that rota- 
tion and hydrodynamical processes do not affect neutrino fluxes either. Indeed, in the region 
where statistically significant deviations from the SSM predictions are seen, helioseismology 
is a much more sensitive probe of the solar structure than neutrino experiments. 

Lithium depletion is another problem for SSM's. The observed surface abundance of 
lithium in the sun (and other main sequence stars) is much lower than the one predicted by 
SSM's. This problem can be solved by anomalous diffusion, which drives lithium inside the 
sun, where it burns. This diffusion could be induced by gravity waves |T^. 



Gravity-wave induced diffusion is a mechanism of mixing inside the solar core. Another 
mixing mechanism, periodical instabilities, has been considered again and again in the lit- 



erature starting with the pioneering work by Dilke and Gough |T6|. Mixing in the solar core 
affects the solar neutrino fluxes, but it should be shown to be consistent with the other data; 
in general, large mixing contradicts seismic observations. Indeed, mixing episodes transport 
from the periphery to the center of the sun not only the desirable ^He or lithium but also 
large masses of hydrogen. The consequent reduction of the molecular weight increases the 



sound speed. This is the main problem of the Cummings-Haxton assumption ^Tj about 



mixing in the solar core: the resulting sound speed is too high [TE\. Brun et al. [^ re- 
cently found that seismic data do not support even small mixing in the solar core. Using 
the turbulent diffusion coefficient suggested by Morel and Schatzman |T^, they obtained a 
sound-speed profile in contradiction with observations. 

The essence of the recent attempts to reconcile the observed neutrino fiuxes with solar 



models pO| , p7| , p!5| , pT| consists in bringing ^He in the solar core by one mechanism or another. 
Arguments against these mechanisms use the accompanying processes, such as the simul- 
taneous bringing-in of hydrogen. The strategy of our work is different: we allow arbitrary 
mixing of ^He, not accompanied by other elements. Moreover, we do not assume that the 
arriving ^He is in nuclear equilibrium with the other elements; an example can be the out-of- 
equilibrium concentration of '^He in mixing episodes caused by instability. Nonequilibrium 
burning of ^He lowers the temperature. We shall allow independent variation of temperature 
within ±5% of its SSM value, assuming that some miraculous fine-tuning keeps the sound 
speed in agreement with seismic observations. In the solar core we also allow radial profiles 
of temperature and ^He density different from the SSM ones. In addition, we vary the two 
relevant cross-sections, 5*17 and 5*34, within very generous ranges. 

We shall demonstrate that the resulting predictions are, nonetheless, incompatible with 
the observed neutrino fluxes. 



II. ASSUMPTIONS, CONSTRAINTS AND NOTATION 

We shall study the production of boron (^B) and beryllium (^Be) neutrinos. Our ap- 
proach consists in comparing the calculated fluxes ($B) ^Be) with the observed ones. 

Inspired by the SSM, we consider an effective Neutrino Production Zone (NPZ) for 
boron and beryllium neutrinos. In the SSM the production of B and Be neutrinos has its 
maximum at r = Rnpz ~ O.O5i?0. When non-SSM profiles are considered, the position of 
this maximum changes. In our calculations we use Rnpz only as a reference radius and take 
into account that the peak of the B-neutrino production does not coincide with that for 
Be-neutrinos. All quantities evaluated at r = Rnpz are denoted by the index NPZ, e.g., 
Tnpz (temperature), pnpz (density), Y^^npz (^He abundance), etc. Temperature, density 
and abundance profiles will be specified. 

We assume that ^He is transported to NPZ by some unspecified process and characterize 
it by the arbitrary value Y^^npz- We do not assume local nuclear equilibrium for ^He: during 
mixing episodes, "fresh" ^He can be brought inside from the outer shells and burned in a 
nonequilibrium regime. However, it is easy to see that thermal equilibrium is established 
very fast and thus ^He has the thermal distribution of the surrounding gas. 

We keep Y^^npz and Tnpz as independent parameters. In a realistic solar model in- 
creasing Y3NPZ increases the ^He + ^He reaction rate more than the rate of the ^He -|- "^He 
reaction, and the temperature must be lowered to keep the luminosity constant; the net 
result would be the decreasing of $Be and $b- In fact, in the SSM one finds that Y3 ~ T~^. 
For our purpose, any dependence between Y^^npz and Tnpz imposes additional restrictions 
and makes our conclusions stronger. 

We do not impose the luminosity constraint. The luminosity sum rule is only used to 
express the pp neutrino contribution to the gallium experiments in terms of the other fluxes. 

The CNO neutrino flux is conservatively neglected to avoid solar- mo del- dependent con- 
sideration. Including the CNO neutrino production would only strengthen our conclusions. 

Helioseismic constraints are used very conservatively. We assume that temperature in 
the solar core can differ from that of the SSM by ±5%, to be compared with a fraction of 



percent allowed for the sound speed squared cj. More specifically, we assume the following 
maximum range of temperature variation: 

0.95 <r^pz/r^fi^< 1.05, (1) 

and the same order of variation for the radial profile of the temperature in the NPZ. 

The sound speed squared is given by c^ = jRT/fi, where 7 and R are constant, T is the 
temperature, and fi the mean molecular weight, l//i = 2X + 3F/4 + Z/2 ^ (5/4) (X + 3/5). 
Therefore, since the accuracy of seismic determination of the sound speed is better than 1%, 
the following constraint must be met: 



X{r) + 3/5 T^SM^^ 



(2) 



XSSM^r) + 3/5 T{r) ' 

Since the allowed variation of T{r) is only a few percent, one can use 

X(r) = X^^*'(r) (3) 

with the accuracy needed: a few percent change in the temperature strongly affects neutrino 
production, while the same variation of X does not do much. 

For the ^He abundance we use Y ^ 1 — X. 

We can use the SSM density profile, p(r), since this profile has been measured by seismic 
observations. 

The astrophysical factors that affect B and Be neutrino production are Su and S34. 

For the reaction p + ^Be the INT Collaboration |Q (1998) suggests the la range 5*17 = 
19+^ eV b; Castellani et al. [||] (1997) suggest Sn = 22.4 ±2.1 eV b. The BBP98 SSM ^ 



and the BTCM SSM JTrf use 19 eV b, while the DS96 SSM uses 17 eV b. We shall use 
the generous range: 

13 < 5i7 < 31 eV b. (4a) 



For the reaction He + He the INT Collaboration [^] suggests the la range 6*34 



0.53±0.05 keV b; Castellani et al. |2| suggest ^34 = 0.48±0.02 keV b. The BBP98 SSM ^ 
and the BTCM SSM Q use 0.53 keV b, while the DS96 SSM uses 0.45 keV b. We 
consider the range: 

0.38 < ^34 < 0.68 keV b. (4b) 

The detected neutrino fluxes used in our analysis are the following. 

For the Chlorine signal we use the Homestake data 2.56 ± 0.21 SNU [^; for B-neutrino 
flux the Superkamiokande data $b = (2.46 ±0.09) x 10^ cm~^s~^ [^; and for gallium signal 
we use the weighted average of the GALLEX || and SAGE @] data 72.5 ±5.7 SNU. 

III. PARAMETERIZATION OF NEUTRINO PRODUCTION RATES 

The rate of Be neutrino production as function of the radius is: 



Qy^Xr) = 47rr2ne(r) m{r) KiiT^r)) (5a) 

= 47rr2 Pj,,{r) n^{r) n^{r) \^^{T{r)) , (5b) 

where n^ {x = 1, 3, 4, 7, e) is the number density of particles (protons, ^He, ^He, ^Be nuclei, 
electrons), Pbc = K7ne/[\e7ne + Xnni] is the electron capture probability and Xij = {vaij) 
is the reaction rate averaged over the Maxwellian distribution of the relative velocities, which 
is a know function of temperature and does not depend on solar model. 

Since the ^Be electron-capture-lifetime in the solar core is very short, TeiJ Be) ~ 1 yr, 
the equilibrium value nj = n^n^Xu/ {neXej + niXn) has been substituted into Eq. ( pa) ) to 
obtain Eq. (|5BD . With Pee ~ 1, the rate of Be neutrino production becomes: 

Q.s.ir) = ^ (1 - ^(0) [rpir)f A34(T(r)) Y,ir) . (5c) 

After similar transformations, the rate of B neutrino production, 

Q.,{r) = Anr'mir) n,{r) X,,{T{r)) = Q..e ^^^4444^^' (6a) 

neyr) Ae7{I [r ) ) 

becomes 

._Nl 2il-Xir))Xir) , XuiTir))X,r{T{r)) 

^us[r) 3 1+X(r) ^^^^J Ae7(T(r)) '^ ^ ' ^^^^ 

Let us introduce the following scaling variables: 

X = t/Rnpz (7a) 

p{x)=p{r)/pNPz (7b) 

y^{x) = {l-X{r))/ (l-X^Pz) (7c) 

y3{x)=Ysir)/Ys,NPZ (7d) 

A34(x) = XuiT{r))/ Xu{Tnpz) • (7e) 

The reference radius Rnpz is taken within the Neutrino Production Zone (NPZ); in practice, 
we use Rnpz = 0.0535/?©. 

The profile function normalizations are p(l) = 1/4(1) = ^3(1) = A34(l) = 1. The total 
rate of Be neutrino production becomes: 

Qi^Be = —^ Rnpz Pnpz{^ - ^npz) nssM ^m{Tnpz) Yz.npz Ipe , (8) 

where we have defined the dimensionless integral over the profile functions 

lBe= dx[xp{x)] y4{x) Au{x) y3{x) , (9) 

Jo 

and A34 is calculated with the SSM temperature profile (the dependence on the cross section 
has been explicitly taken into account by 5*34). 

In practice, we can use Pnpz = pffpz ^^^ ^npz = -^jvpz; and in the integral 
p{x) = p^^^{x) and 2/4(0;) ~ yf^'^'^ix), since these quantities are strongly constrained by 



helioseismology and the dependence of the results on their precise value is not strong. On 
the contrary, A34(a;) and 1/3 (x) are not directly constrained by helioseismology and can be 
different from the SSM: therefore, we shall study the effects of their change. Therefore, the 
Be neutrino production rate can be expressed as 

QuB^ _ 'S'34 ^u{TnPz) y^,NPZ Ibc /^P n 



nSSM qSSM \ frpSSM\ vSSM tSSM ' 

^UBe '-'34 '^^^K-'-NPZ) ^3,NPZ ^Be 

Notice that A34(T) is a known function of T, independent of the solar model. 
Similarly, the B neutrino production rate can be expressed as 

QuB _ 'S'34 5*17 rjiTNPz) y'i,NPZ Ip (^(\U'\ 



nSSM qSSM qSSM „(rpSSM\ \^SSM tSSM 

^UB '^34 "^l? Vi-'-NPZ) ^3,NPZ -', 



!uB '-'34: ^17 'IK-^NPZJ ^3,NPZ ^B 

where 

,(.,.MgMI) ,^ 

is another known solar-model-independent function of T, and 

Ib= dx[xp[x)\ — H{x)y3{x) 12) 

Jo 2 — y4:{x) 

where H{x) is another profile function: 

viTir)) ^ viTNPz)H{x) , (13) 

whose normalization is also H{1) = 1. 

At this point, we should discuss the role of the normalization radius Rnpz- Equa- 
tions (|T0|) show that the only dependence on Rnpz left is in the profile functions /. As 
long as we compare models with the same profiles, the functions / cancel and there remains 
no dependence on Rnpz- However, we are interested also in models with non-SSM pro- 
files. In this case variations of the reference radius modify, in principle, the ratios I j I^^^ 
and consequently the fluxes, when all other parameters stay fixed. The important point is 
that these changes of fluxes and ratios can be absorbed by appropriate variations of other 
parameters, in particular of Tnpz and Y^^npz, and therefore, a reference radius does not ap- 
pear as another free parameter. In particular, we use such a parameterization of the profile 
deformation (see below) that a variation of the reference radius is equivalent to an overall 
multiplicative factor for the quantity whose profile deformation is considered. 

It is also clear that the choice of the SSM is not essential: Eqs. ([T0|) are valid for any 
two solar models that satisfy the helioseismic constraints. 

1. SSM radial profiles 

We start our analysis by assuming that the shapes of all radial profiles, including T{r) 
and Y^{r), are those of the SSM. Since the SSM gives a good description of the solar structure 

6 



and since we rescale all values at Rnpz, this very reasonable assumption should already be 
a good approximation. In this case I Be = -^Jf^^ and Ib = Ib^^^ in Eqs. (10). Then one 



obtains 



17 



QvBe. 'S'34 / Tnpz \ Yz^NPZ 



X ^Tssir (14a) 



nSSM qSSM ^^ rpSSM I ^^ vSSM 

^UBe "-^34 X-'-NPZ / ^3,NPZ 

QuB S34 ^^ Sn [Tnpz\ y3,NPZ /-i^, X 



nSSM qSSM qSSM rpSSM I vSSM 

^ub '-'34: ^17 K-'-NPZ / ^3,NPZ 

where, for convenience and without loss of generality, the (model-independent) temperature 
dependence of the two functions A34(T) and ^(T) have been approximated by power-law 
functions (this approximation is quite good in the region of interest T^pz/T^pz = 1±0.05): 
A34(T) ~ T^^ and ?7(T) ~ T'^". We shall use these power-law functions from now on. 

At this point we are left with four free parameters 5*34, S'i7,T/vpz and Y^^jypz- Eqs. (|T3p 
show that changing Y^^^pz by an overall factor has the same effect as changing 6*34. Since 
Up and UBe fluxes are proportional to Y^^pz, increasing ^3 tvpz results in the simultaneous 
increase of both fluxes without changing the ^pe/^B ratio. It does not help to solve SNP, 
making actually the problem even harder. 

In Figure |I| the two thin solid curves on both sides of the curve "temp." (thick solid 
curve) bound the allowed region. This region is given by Eqs. (|l^) when Sn and 6*34 are 
varied within the limits of Eqs. (^ and the temperature is arbitrary. In this region the best 
fit to the experimental data is indicated by the small cross labeled 8 ($Be = 2.02 x 10^ 
cm-2s-\ $B = 2.34 X 10^ cm'^s'^) and it is obtained with ^34 = 0.38 keV b, ^17 = 31 eV b 
and Tp^pz/T^p^ = 0.969. Even if we are assuming no CNO flux, its x^ = 35 is still quite 
large: more than 5a away from the experimentally allowed region. It is easy to understand 
from Fig. [| that further lowering the temperature below T^pz/T^pz = 0.969 does not help, 
unless one also allows values of Sir larger than 31 eV b (larger upward shifts). 

It is instructive to follow some selected trajectories within the allowed region. 

We can start for example from the BBP98 SSM, i.e., the diamond in Figure |I]. 

When the temperature is lowered and the other parameters are kept constant, the fluxes 
follow the thick solid line labeled "temp." [|. The heavy dots labeled 1 and 4 correspond to 
temperatures scaled by 0.98 and 0.95. From these points 1 and 4 one can change 5*34 (or 
equivalently Y^^ssm) and reach points 2 and 6, respectively. Then one reaches points 3 and 
7 by increasing S'17. Trajectories are shown by dashed lines. Going from point 1 to point 
2 corresponds to changing 6*34 from 0.53 to 0.38 keV b, and going from point 2 to point 3 
to changing 5*17 from 19 to 23 eV b. In fact, this set of parameters is the best solution, 
if one keeps 5*17 within the range of Eqs (^) and the temperature within 2% of the SSM 
value. Its x^ = 48: it excludes the experimental data at more than 6a. For reference, had 
we taken ^cno = ^c^no^ instead of ^cno = 0, the x^ would have been 82, which excludes 
the experimental data at more than 8a. 



^This is not the usual temperature solution, where the other parameters, and in particular the 
^He abundance, also change as function of temperature. 



Two other possible trajectories corresponding to the lower temperature (Tnpz/T§p^ = 
0.95) are the one leading to point 5 (x^ = 80), which is reached by increasing Sn from 19 
to its maximum value 31 eV b, and the one leading to point 7 (x^ = 187), which is reached 
by decreasing 5*34 from 0.53 to its minimum value 0.38 keV b (4 — > 6) and then increasing 
5*17 from 19 to 31 eV b. Note, that we described these trajectories only for illustration; 
numerical calculations used directly the parameterizations of Eqs. (0). 

2. The modified radial profiles 

Now we generalize the above analysis and allow the shape of the radial profiles to be 
different from that of the SSM. Therefore, we do not assume anymore lBe,B = ^bc'e ^^"^ 
obtain: 

r) /nSSM cSSM / rp \ -13 J ijSSM 

In comparison with the previous case the only remaining hope is that the ratios 
I Be,B / Ib^^b 1 determined by the radial profiles, could improve the agreement with the ex- 
perimental data. 

In particular, a large deformation of the ^He profile could, in principle, significantly 
change the integrals lBe,B circumventing the problem of the too large ratio of boron to 
beryllium flux; this approach corresponds to the proposal of Gumming and Haxton |T7|. 



In fact, the only profile functions which can significantly affect the integrals I Be and 
Ib are those of the temperature and the ^He abundance: the former because of the strong 
temperature dependence of the rates, the latter because of its being basically unconstrained 
by helioseismology. 

To estimate the dependence of these integrals on the temperature and ^He profiles we 
introduce the following parameterization: 

FMR,,.) - ^sMrla...) X y^£^^^f . (16) 

where F is either the temperature or the '^He profile. This parameterization keeps the nor- 
malization -F(l) = 1 fixed, gives larger deformations for larger |(5|'s and reproduces the SSM 
profile when 5 = 0. Positive (negative) (5's increase (decrease) the function for r > O.O535i?0 
relative to r < O.O535i?0. Since we want a single parameter function, the scale O.O535i?0 
is kept fixed: 5 is sufficient to control the first derivative of the profile at Rnpz- In other 
words, this parameterization includes all possible changes of the overall normalization and 
first derivative of the temperature and ^He profile. Parameterizing differently the deforma- 
tion in the NPZ (more parameters and/or different functional forms) must lead to similar 
conclusions, since the beryllium and boron neutrino production region is sufficiently small, 
so that higher order derivatives result only in small corrections. 

In fact, we have explicitly checked that, if we drastically change this scale from 0.05Rq 
to O.O25i?0 (changes in the other direction are not relevant to the solution of the solar 
neutrino problem), this new parameterization reproduces the same results within a few 
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percent, obviously for different values of (5's (for instance, the case S = —4 is reproduced 
hy S = —2.7). The dotted curve in Fig. |^ shows that the difference between the two 
parameterizations is small in the region where neutrinos are produced. In this respect, we 
find instructive, and less dependent on the specific parameterization, to consider the changes 
of the profile at two representative points r = O.Oli?© and r = O.IRq of the inner and outer 
part of the NPZ instead of 6 itself. 

Notice that changing the reference radius Rnpz — > R'npz i^ Eq. (|TB|) only give an overall 
multiplicative factor that can be absorbed in the scale factors T/vpz and Ys^atpz, confirming 
that the choice of Rnpz is irrelevant also in the case of deformed profiles. 

Temperature profile with \5\ = 0.057 in Eq. (|T6|) results in a maximal temperature cor- 
rection of about 5% in the relevant region {O.OIRq < r < O.IRq). In Table | we report the 
dependence of Ibs and Ib on the change of shape of the temperature profile. It is particu- 
larly important the relative change of Ibc/Ib- The ratio Ibc/Ib changes at most of about 
15%, relative to the SSM ratio. Changes of this magnitude of the i>b and ube fluxes are 
compatible with their theoretical uncertainties (about 10% for the z/^e and about 30% for 
the ub flux). Therefore, modifying the temperature profile does not help much. 

The ^He profile affects the neutrino fluxes much more strongly, because the ^He abun- 
dance 1^3 (r) is not constrained by helioseismology and is allowed to be considerably different 
from the SSM shape. One can consider in this case larger values of \6\. 

In Table |I| we report the dependence of Ise and Ib on the ^He profile. Note that 
increasing the abundance in the center relative to the outer part of the NPZ {6 < 0) has 
the effect of decreasing the ratio of the Be to B neutrino flux. It is possible to suppress this 
ratio by a factor almost 0.6 {6 = —4) at the price, however, of increasing (decreasing) the 
inner (outer) part of the NPZ by a factor about 8 (1/4). Such a change of profile (a factor 
of thirty over less than O.IRq) is really dramatic in comparison with the SSM as one can 
see in Fig. 0. 

In the following we analyse why even large deformation of the ^He profile can only 
partially improve the comparison with the experiments. A deformation of the ^He radial 
profile that increases the ^He abundance for r < O.O535i?0 and decreases it for r > 0.0535Rq 
boosts $B relative to ^Be, since the boron neutrinos are produces at smaller radii than 
beryllium ones. However, such a deformation has the effect of moving both production 
regions to smaller radii with two consequences: both fluxes increase (even if by different 
factors) and the difference between the two production regions shrinks. On the one hand, 
the more we deform the profile the more is difficult to keep boosting ^b relative to ^Be'- this 
explains why even the large deformation we have considered {6 = —4 in our parameterization, 
which corresponds to an increase of ^He in the inner part of the NPZ relative to the outer 
part by more than one order of magnitude) can only improve the ratio by a factor 0.64 
(see Table 0). On the other hand, even if a strong deformation results in the required 
small I'Be/^B ratio, it would be still necessary to strongly reduce both fluxes. This could 
be achieve either by lowering the temperature or by lowering the overall ^He abundance 
(the normalization factor Y^^pz)- The possibility of lowering the temperature is limited 
by helioseismology, and it does not seem physically possible to strongly reduce the ^He 
abundance in the solar core (less strongly towards the center so that the required profile 
deformation is achieved), since any mixing should increase the ^He in the core. 

In Fig. 1^ the thick solid line shows the variation of neutrino fiuxes with temperature in 
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the case of the strongest deformation of ^He profile {S = —4). The thin sohd hnes bound 
the region that is spanned when also Sn and 5*34 are varied within the limits of Eqs. @. 
The best fit point is labeled 4 {^Be = 1-72 x 10^ cm-^s"^ and $b = 2.41 x 10^ cm-^s-^); 
it is obtained with 5*34 = 0.38 keV b, ^17 = 31 eV b and Tnpz/TNPZ^^^ = 0.95 and has 
a x^ = 32, or about 5cr's away from the area allowed by experimental data. The trajectory 
2 ^ 3 ^ 4 shows a possible way to reach the best-fit point 4. 

Including contributions of CNO neutrinos to the signals makes the agreement with ex- 
perimental data even worse. 

We could not reach in our calculations the point of Gumming and Haxton (shown in Fig. |I| 
and Fig. ^ by an asterisk) because of the seismic constraints imposed in our calculations 
and because we do not allow an overall strong reduction of ^He in the core, which cannot 
be caused by mixing. 

IV. CONCLUSIONS 

The solar neutrino experiments have the status of disappearance oscillation experiments. 
This statement is based on the impossibility of explaining the observed deficit of neutrino 
fluxes with astrophysical processes not included in SSM's and revised nuclear cross-sections. 
Helioseismic data strongly constrain the possible non-SSM astrophysical processes. 

Recently the idea of ^He mixing in the solar core has been revived ||l5| , |20| , pl| , [T7|| . The 



abundance of ^He is not directly constrained by seismic data. In principle this fact opens 
a road to possible revisions of the SSM predictions for the neutrino fluxes. However, in 
any realistic model, ^He mixing is accompanied by other phenomena, e.g., hydrogen is also 
transported into the core; because of these accompanying phenomena ^He mixing can be 
constrained by seismic observations (see |T8| , p!T| ). We present a more general approach where 



constraints on the neutrino fluxes are valid for any mechanism of '^He mixing. 

We assume arbitrary ^He mixing. Some unspecified process brings fresh ^He from the 
^He-richer outer shells into the solar core where thermal equilibrium is quickly established, 
but not nuclear equilibrium. Since this process could consist of short mixing episodes, we do 
not impose the solar-luminosity constraint. The density radial profile is taken from seismic 
data. Within the accuracy needed for calculations of neutrino fluxes, the X{r) and Y{r) 
profiles are also provided by seismic data. Therefore, the only solar parameters that are 
needed for model-independent calculations of neutrino fiuxes are the temperature T, the 
^He abundance F3, the astrophysical factors S'34 and S'17 (cross-sections), and the radial 
profiles of T{r) and 1^3 (r). We allow independent and large variations of T, S'34 and 5*17 as 
given by Eqs. (|I|) and (^. 

Regarding the radial dependence of the temperature, T{r), and of the ^He abundance, 
Y^i^r), we used two approaches. 

The first one consists in using the SSM radial dependences. This choice appears reason- 
able, because B and Be neutrinos are actually produced in a narrow region and distortions 
of the radial profiles, if not extreme, should not change the fiuxes much. Results of this first 
approach are presented in Fig. |I[ 

The thick solid line describes the evolution of the fiuxes with temperature. The two 
thin solid lines confine the region allowed by arbitrary variations of '^He abundance and the 
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temperature, accompanied by variations of 5*34 and 5*17 within the ranges given by Eqs. (^. 
The cross labeled 8 shows the best fit, which has a very large x^ = 35. 

In the second approach we also allow changes in the shapes of the radial profiles T{r) and 
Kj(r). We find that the distortion of the ^He profile has a much stronger effect on the neutrino 
fluxes and it is the only important one in our analysis. To improve the agreement with 
experimental data, the ^He abundance should, contrary to the SSM case, increase towards 
the center, as illustrated by the dashed line in Fig. |^. We are aware of no physical mechanism 
in solar models that could produce such a dependence. We just take an unconvential and 
very strong radial dependence such as the one shown in Fig. ^ as an ad hoc assumption. 

With this extremely deformed ^He-density profile we repeat the exercise of varying the 
^He abundance and temperature in the neutrino production zone, with simultaneous varia- 
tions of Su and Sn. In Fig. | the allowed region is confined by the two thin solid curves. 
The best fit to observational data is given by point 4, and it still corresponds to a large 
X^ = 32, i.e. more than 5a away from the experimentally allowed region. 

In all our calculations we have conservatively assumed a vanishing CNO neutrino flux. 
Inclusion of any amount of CNO neutrinos makes the disagreement with the experimental 
data worse. 

We conclude that practically arbitrary ^He mixing, which also includes physically unjus- 
tified distortions of the ^He radial profile, and independent variations of temperature and 
cross-sections cannot explain the observed solar neutrino fluxes. 

One might ask, however, whether and to what extent the considered uncertainties, espe- 
cially in ^He mixing, can reduce the discrepancy between the SSM solar-neutrino fluxes and 
observations. The most important source of uncertainty is the p + Be cross section (S'17) , 
which affects only the predicted flux of B neutrinos. Taking S'17 40% lower than presently 
used and the core temperature T^ 1.4% lower (as maximally allowed [^ by seismic obser- 



vations) one arrives at the minimum B-neutrino flux 3 ■ 10^ cm^^s^^ |2^, i.e. 7.4cr higher 
than the measured one. It is more difficult to estimate a reasonable effect of ^He mixing 
(if this process exists at all). It is constrained by accompaning processes, such as hydrogen 
mixing, and a realistic model is needed for such calculations. The method used in this work 
can give only an upper limit of the influence of ^He mixing on neutrino fluxes. The ad hoc 
assumptions used in this analysis are maximally favorable for this influence, though rather 
unrealistic. 
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FIGURES 

FIG. 1. Neutrino fluxes allowed by arbitrary ^He mixing accompanied by independent varia- 
tions of temperature, 834 and Sn (the ^He and temperature radial profiles are those of the SSM's). 
The horizontal (vertical) axis shows the berylhum (boron) flux; fluxes are measured in units of the 
reference SSM of Bahcall and Pinsonneault BP95 [^ on the bottom (left) scale, and in cm~^s~^ 
on the top (right) scale. The region allowed by the SSM's (BP95 SSM 90% confidence region) is 



shown by the dotted ellipse; the diamond shows the BBP98 SSM |10| prediction with the relative 
la errors and the square shows the DS96 [^] SSM. The solid ellipse confines the region allowed at 
3(7 by the experimental data (if there is any contribution to the signals from CNO neutrinos this 
region becomes smaller). A representative set of nonstandard solar model calculations are shown 
by small "x" symbols and is taken from the review by Hata and Langacker [^, while the asterisk 
indicates the ^He mixing assumption by Cunning and Haxton |17|. The thick solid curve "temp." 



gives the evolution of the predictions of the BBP98 [[10[| model with variations of temperature 



according to Eqs. (|14|). In particular points 1 and 4 correspond to temperatures lower than Tssm 
by factors 0.98 and 0.95, respectively. When Sn diminishes, points move vertically down; when 
5*34 diminishes, they move diagonally towards the origin. Three trajectories (1 — > 2 ^ 3, 4 ^ 5 
and 4 — 5- 6 — > 7) are shown for illustration. The two thin solid curves bound the region allowed 
by arbitrary changes of the '^He density and of the temperature, accompanied by variations of ^34 
and 5i7 within the limits given by Eqs. (HI). The best fit to the experimental data is labeled "8" 
and has the (too large) x^ = 35. 

FIG. 2. Profiles of the ^He abundance normalized to 1 at Rnpz = 0.0535Rq. The solid curve 
shows the profile in the SSM, while the dashed one shows the profile corresponding to the largest 
deformation considered in this paper, i.e., 6 = —4 in Eq. (p^) and in Table |II|; the dotted curve is 
an example of different parameterization as discussed in the text. 

FIG. 3. Same as Fig. |||, but with a modified radial profile of the '^He abundance. This case 
corresponds to the maximal distortion {5 = —4) of the SSM profile considered in this work (see 
Table || and Fig. ^. The temperature trajectory (towards lower temperatures) is shown by the 
thick solid curve "temp.". The points 1 and 2 correspond to temperatures lower than Tssm by 
factors 0.98 and 0.95, respectively. The trajectory 2 ^ 3 — s- 4 is shown for illustration: along the 
track 2—5-3 ^34 decreases by a factor 0.38/0.53, while along the track 3^4 5i7 increases by 
a factor 31/19. The two thin solid curves bound the allowed region. For comparison the region 
allowed when keeping the SSM radial profile (see Fig. [l|) is bound by the two dotted lines (thin 
solid lines in Fig. |^). The best fit is given by point 4 {x^ = 32), more than 5a away from the 
experimentally allowed region. 
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TABLES 



TABLE L Dependence of the integrals Ibc and Ib, Eqs. (^) and (12), on the temperature profile, 
parameterized according to Eq. (p^). The temperature profile is deformed keeping T{Rnpz) fixed. 
The first column shows the parameter 6 used in Eq. (|l6|). The second and third columns give 
the corresponding change of the integrals I Be and Ib, while the fourth column shows the change 
of their ratio. The last two columns illustrate the tilting of the profile relative to the SSM by 
reporting T/T^^^'^ at r = O.Oli?© and r = Q.IRq. 



IbJI'bI'' 



Ib/I'b'^' 



{lBe/lB)/{lBe/lB) 



SSM 



T/T^^^'{0.01Rq) T/T^^^\<d.lRc^ 



-0.057 


0.95 


1.21 


0.79 


-0.044 


0.96 


1.14 


0.84 





1.0 


1.0 


1.0 


0.044 


1.08 


0.94 


1.15 


0.057 


1.11 


0.94 


1.18 



1.03 


0.98 


1.02 


0.98 


1.0 


1.0 


0.98 


1.02 


0.97 


1.02 



TABLE II. Dependence of the integrals I Be and Ib, Eqs. (^) and (0), on the profile of the 
^He abundance. The information is analogous to the one in Table |. 



T I tSSM 
^Be/^Be 



IB/If'' 



{lBe/lB)/{lBe/lB) 



SSM 



ys/yl'^'m 


]li?0) 


ys/yf''^{0.lRe 


8.06 






0.24 


4.78 






0.34 


2.84 






0.49 


1.69 






0.70 


1.0 
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0.59 
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0.12 
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1.17 


2 


1.32 


0.93 


1.41 


3 


1.66 


0.96 


1.72 


4 


2.23 


1.03 


2.16 



15 



CO "n 

(D 1-' 



H%/H%SSM 



o 
b 



o 



3 



dd 

CD 



dd 

CD 
w 









CO \ / 




# 


m 


♦ 


X 






l\ 




"^ / 
















\\ 




S 17 




o 


a 


CO 


g 


- 




l\\ 






el 


M 


tn 


o 






"\\ \ 




X ^ 




3 
3 


CO 


13 


a 






\\\ 




V 




0) 


CD 
00 


CO 


■ 




\ \ \ 








5 


M 


M 


o 






w \ 








1 


s 










\ \ \ 








w 


"n 





a 






■ \\\ 




# 




5» 
i< 
rt- 
O 
P 


i 


o 




- 




\ \ ''^ 


\:^ 






3 

b: 






^ 






W 


A 
















\ \ 




\ 











S 






\ V 




\, 




OQ 






a 








\- 




M 






sr 






\ \ 


4^ 


ai\. CO 




" 






3 






\ 


K 






5' 







O 


■ 


o 


\k 




\ ^ \ 










ii 




en 




X 


\ GO \ 
\ X 

V '* 

\ ** ^\ 

\ '^^V 


1 

• \ 




X 


X 


- 


- 






3 


\ 


X 


r 






- 






£- 


\ " 
























■ ^^\ * 


































\ 


1 








1^ 


■ 




1 








1 C/2 








\ 










S 


)( 




\ 






1 1 


^ 


, , , , \ , 


1 






— 1 — 


1 





CO 



hf^ 



8 



(|)(^B)[10Vm-^s-^] 



05 



dd 

^CD 



o 

CO 

O 



00 M 

CO 



- 4^ 



Oi 



16 



CO T1 
<D 1-' 
-) -Q 
<D ■ 
tl IV) 



He abundance 



3 



o 




17 



CO "n 

(D 1-' 



CO 



H%/H%SSM 



o 
b 



o 



3 



dd 

CD 



dd 

CD 
w 





I 








CO \ / 


1 




# 


■ 


♦ 


X 






1\ 








"^ / 


















\ \ 








S 17 






o 


a 


CO 


g 


- 




\ '^ 












el 


M 


tn 


o 






'\ ■ 








X ^ 






3 
3 


CO 


13 


a 






\ 








V 






0) 


CD 
00 


CO 


■ 




A 














5 


M 


M 


o 






A 














1 


s 










' \ 


\ 












w 


'to 





a 








\\ 


\ 


\ 


# 






5» 
i< 
rt- 
O 
P 

3 

b: 


i 


o 




■ 








\ 




03 \, ^ 






s 






a. 


- 










\ ^- - ^^ 






M 






^ 








X 




\ 


^^ 






13 






a 












\ 


X ^ 


\ 








P 
O 

re 


- 












x\. 


^^^ 




5' 






1 


- 












\ w 


^^ 




B 




















^a X 












10 










. X 


\ 






^^^ 


h-* 






s 


- 


o 








\ 


J^v 




^ 


^ J2 










en 








x\ 

3 


\ 

\ X 

\ X 

X '. 1^ 


X 

\ ^ 

a 

X 

X 

X \ 


X 

\ 

1 


X 
X 


*< 


X 


1^ 


N 




1 








1 C/2 
















\ 








S 


)( 












— 1 




^ 


1 




\ 






— 1 — 


1 





CO 



h^ 



dd 

^CD 



o 

CO 

o 



00 M 
CO 



- 4^ 



Oi 



05 



8 



(|)(^B)[10Vm-^s-^] 



